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ABSTRACT 



Capacitance and conductance measurements using two and 
four wire techniques were developed and statically and 
dynamically calibrated in this thesis. The voltage 
sensitivities range from 7.3 to 8.1 ± 0.1 mV/cm for the two 
wire capacitance system static calibrations. This is ± 5.2% 
of the limiting theoretical value. The voltage sensitivities 
range from 0.3 to 0.4 ± 0.1 V/cm for the four wire 
conductance system static calibrations. Dynamic calibrations 
were only completed for the conductance system. The dynamic 
calibration results were weakly frequency dependent with a 
co-0.15 decay in a limited, 2-4 Hz range. 

Wind power spectrum measurements were taken in the 
existing Upper Ocean Simulations Facility at the Naval 
Postgraduate School . There was excellent agreement in the 
spectra with both techniques. 

Driven gravity wave frequency downshifting and wind 
energy dumping was observed in the combined gravity wave and 
wind-wave measurements. The power spectra peaked near two 
Hertz and decayed at 50 to 70 dB per decade, or as co“5 to 
co-7 f or both systems. 

Gravity wave phase speed and wavelength measurements 
were performed with the capacitance system. The results were 
approximately 40% higher than theory. 
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I . INTRODUCTION 



A. BACKGROUND AND MOTIVATION 

In the 1980's, there were tremendous advances in the 
theories that describe nonlinear systems in many fields of 
science and engineering [Ref. 1]. These nonlinear systems 
vary in complexity from a simple driven pendulum, to global 
weather systems. This thesis studies the nonlinear system of 
wind wave interaction on a "deep ocean" using an existing 
wind wave tank, Upper Ocean Simulations Facility at the 
Naval Postgraduate School (NPS) . 

A theory describing the power spectrum of ocean surface 
wave crests, independent of the fetch (the distance over 
which wind blows over water) and duration of the generating 
wind, was developed by Phillips [Ref. 2]. The Phillips 
spectrum, based on dimensional analysis, shows a © - 5 
spectrum for large angular frequencies, co . 

Zakharov and Filonenko later described the stochastic 
oscillations of the surface of a liquid with a © - ^ power 
law spectrum [Ref. 3]. In the late 1980's, Larraza applied a 
© - 4 power law spectrum to wind waves, flicker noise, solar 
wind, and classical second sound [Ref. 4]. 

There is a thirty year history of reported field and 
tank wind wave and driven gravity wave spectrum 
measurements, summarized nicely by Donelan, et . al [Ref. 5]. 
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When Dr. Larraza came to the Naval Postgraduate School, we 
had an opportunity for experimentalist and theorist to 
collaborate in the laboratory while exploring his recent 
power law spectrum and second sound theories in a controlled 
environment . 

There is also a history of laboratory controlled, in 
wind wave tanks, wind driven phase speed measurements. Some 
recent experiments which have been done to understand this 
phenomenon has been either flawed [Ref. 6] by aliasing due 
to excessive spacing of probes or were not designed to look 
for the effects of a turbulent restoring force [Ref. 7, 8]. 

Wave turbulence are those processes in open driven 
systems which are dominated by a random redistribution of 
energy (generally due to time reversible nonlinear 
processes) with effectively infinite degrees of freedom 
[Ref. 9]. Significant theoretical progress in the study of 
random wave systems driven far off equilibrium has lead 
Larraza, Garrett, and Putterman to the startling prediction 
of a second surface wave mode [Ref. 10]. This "second sound" 
mode corresponds to propagating compressions and refractions 
of the surface roughness. 

B. SCOPE 

This is the first thesis involved in the study of 
nonlinear wave turbulence using the existing Upper Ocean 
Simulations Facility at the Naval Postgraduate School. The 
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primary purpose of this thesis is to develop, test, and 
calibrate a wave height measurement system that may be used 
in subsequent research efforts studying various wave 
turbulent effects. 

There were two wave height measurement system 
investigated in this thesis. A two-wire capacitance system 
and a four-wire conductance system were independently 
developed and tested. 

Preliminary measurements were made to verify the 
dispersion relation for driven deep water gravity waves 
[Ref. 11] 

0)2 = gk, 

Eqn . 1.1 

where g is the acceleration of gravity and k is the wave 
number, and an investigation of the wind wave power spectral 
density and quality factor Q of the tank were conducted to 
help identify any tank limitations. As it stands, the second 
sound theory requires a power spectrum that decays as co - ~*. 
All these measurements are designed to lead to successive 
experiments which may confidently measure wind wave second 
sound . 

This thesis has four chapters which describe the various 
aspects of the experiment. Chapter II discusses the theories 
of the physical properties which will be measured, and the 
theories of the various measurement techniques used. Chapter 
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Ill describes the wind wave tank facility and it's 
capabilities, the experimental apparatus construction, 
development, and the various calibration techniques for both 
systems, and the various experimental measurement techniques 
used. Chapter IV will present measurement results. Finally, 
Chapter V will discuss the conclusions of this work, and 
suggest directions for future research efforts. 
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II. THEORY 



A. INTRODUCTION 

This chapter discusses theories related to gravity- 
waves, in three broad categories. These categories are 
variables, measurements, and techniques. The chapter will 
define surface gravity waves and energy density variables 
and then derive surface gravity wave dispersion 
relationships. The theories for cylindrical capacitance and 
four wire measurements will then be discussed. Finally, 
least squares fit and deviation data reduction techniques, 
Hewlett-Packard (HP) cross correlation measurement 
techniques, phase speed measurement, wavelength, and quality 
factor measurement techniques will be discussed. 

B. SURFACE GRAVITY WAVE AND ENERGY DENSITY VARIABLES AND 

DISPERSION RELATIONS 

The free surface of a liquid in equilibrium in a 
gravitational field is a plane. If, under the action of 
some external perturbation, the surface is moved from 
its equilibrium position at some point, motion will 
occur in the liquid. This motion will be propagated over 
the whole surface in the form of waves, which are called 
gravity waves, since they are due to the action of the 
gravitational field. Gravity waves appear mainly on the 
surface of the liquid; they affect the interior also, 
but less and less at greater and greater depths. [Ref. 
12 ] 

A reciprocal relation between a deep water surface 
gravity wave phase speed and frequency will be derived in 
this section. Deep water is defined as a semi-infinite fluid 
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below a free surface. The fluid depth H is much greater than 

X 

the reduced wavelength ~ . An intuitive feel for the surface 

271 

gravity wave dispersion relation was developed by Sir M. J. 
Lighthill, using simple energy arguments (Rayleigh's 
method) [Ref. 13] . A more rigorous derivation using invisid 
equations for incompressible flow in a constant 
gravitational field can be found in "Fluid Mechanics" by L. 
D. Landau and E. M. Lifshitz [Ref. 14]. Both developments 
will be presented here. 

1. Rayleigh's Method 

The gravitational potential energy PE per unit 
horizontal area, at any point on the surface, is found by 
integrating the potential energy per unit volume from the 
surface to the bottom. The maximum excess potential energy 
per unit horizontal area over the displacement of the 
undisturbed free surface (^ = 0}is 

PE = \ pgC 2 , 

Eqn . 2.1 

where p is the volume density of the liquid. 

The kinetic energy KE per unit horizontal area of a 
sinusoidal wave is 
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KE = 2 P k ' 1( |^ )2 - 



Eqn . 2.2 



where t is the time 

The mechanical energy E per unit area of the system 
is E = KE + PE. The deep water surface gravity wave 
dispersion relation is found by setting the maximum value 
of equation (2.1) equal to the maximum value of equation 
( 2 . 2 ), 

©2 - gk ; 

Eqn . 1.1 

which can be further reduced using phase speed c, c = ©/k, 
and the relationship between angular frequency and natural 
frequency, © = 2nf , to 




Eqn . 2.3 

The above simple result assumed that there was an 
exponential decay of velocity away from the surface which 
restricts equation (1.1) to deep water. 

2. Invisid Equations 

The derivation using invisid fluid equations begins 
with Newton’s Second Law, F = ma. Bold letters represent 
vectors. The equation of motion of a fluid is determined if 
this force is equated to the gradient of the pressure p on a 
unit volume, F = ma = Vp. The equation of motion of a fluid 
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was developed by L. Euler in 1778 [Ref. 15]. Euler's 
equation in a gravitational field is given by 



5v 

5t 



(v* V) v 



V£ 

P 



+ g. 



Eqn . 2.4 



where v is the particle velocity. 

Fluid particles in a gravity wave travel on the 
order of the amplitude A of the wave during a time interval 
of period T. The particle velocity is of the order A/T. The 
velocity varies noticeably over time intervals of the order 
T and distances of the order of the wavelength X in the 

direction of propagation. The time derivative of the 

5v _ 

velocity ~ is of the order v/T, and the space derivative Vv 



is of the order v/X. The left side of equation (2.4) is 
approximately equivalent therefore to 



5v 

5t 



(v V) v 



A 1 1 fA'te 

T T + X VT/ 



Eqn . 2.5 

We can neglect the (v*V)v term if the amplitude of the 
oscillations in the gravity wave are small compared with the 
wavelength, A << X, and we have potential flow. 

The equation of continuity in a potential flow of an 
incompressional fluid is 

v 2 4) = 0, 

Eqn . 2.6 
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i.e. Laplace's equation for scalar potential <{>, v = V<{>. 
Integrating equation (2.4) after neglecting the (v*V)v term 
and a little algebraic manipulation, gives 



56 

P = -P9Z - P 5t 



Eqn . 2.7 

Let C, = £(x,y,t) be the z coordinate of a point on 
the surface. If a constant pressure p Q acts on the surface, 
equation (2.7) reduces to 



0 



gC + 



Ml 

5t' z =;- 



The z coordinate of a point on the surface is now 



Eqn . 2.8 



c = 



1 6 £| 

9 5t'z=C‘ 



Eqn . 2.9 

The vertical component of the velocity of a point on 
the surface is 




56 

but v z = so that 



M) = K 

&z)z=C > 6t 



Eqn . 2.10 



Eqn. 2.11 
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Substituting equation 2.9 into 2.11 we obtain 




5z 



1 S 2 ^ 

3 5t 2 J 



z=C 



0. 



Eqn . 2.12 

The system of equations which determine the motion 
in a gravitational field are: 



V 2 <j) = 0, 

Eqn . 2.6 



8 & 1 
+ “ 



sV 



5z 



g 8t‘ 



z=C 



0 . 



Eqn . 2.12 

Assume that the velocity decays exponentially away 
from the surface at a rate k. The velocity potential of the 
deep water surface gravity wave propagating in the x- 
direction and uniform in the y-direction is. 



<t> = A Q e^ z cos (kx-cot ) . 

Eqn. 2.13 

Substituting equation 2.13 into the boundary conditions of 
equation 2.12, we obtain 




Eqn . 2.14 
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or 



g) 2 = gk . 

Eqn . 1.1 

This gives the same deep water surface gravity wave 

dispersion relation between wave number and angular 

frequency, as developed using Rayleigh's method. 

The deep water approximation will start breaking 

down when the bottom of the tank interferes with the 

solution given by equation (2.13). This will occur if the 

e term is non-negligible at z = -H, H being the depth of 

1 

the water, which occurs when k < ^. For our tank, H = lm so 



the deep 
_ 2 71 

X » ~j— = 2nH = 



water 



6.3 m and / >> 



approximation 

-S/5 

2n\la 



IS 



= 0.20 Hz. 



good 



for 



C. TRANSDUCER PRINCIPLES 

Understanding the electrical properties that are 
proportional to changes in surface water height are an 
essential part of this thesis. The electrical properties and 
measurement techniques used in this work will be discussed 
next . 

In an electrical circuit, complex electrical impedance 
is proportional to voltage changes when a constant current 
is applied to the circuit. The first type of transducer 
described in this thesis is based on applying a constant 
voltage and measuring the change of capacitance with respect 



11 



to the length of a cylindrical capacitor as a function of 
water height. 

The second type of transducer described in this thesis 
is based on applying a constant current and measuring the 
change in the magnitude of the complex impedance with 
respect to an effective cross sectional area of a volume of 
water under test. The impedance is inversely proportional to 
the water height. 

This section will detail both transduction properties 
and the associated measurement techniques . 

1. Cylindrical Capacitance 

A vertical insulated wire penetrates the surface of 
water approximately forming a cylindrical coaxial capacitor. 
The wire is the inner conductor, the insulation is the 
dielectric, and the water bath acts as an outer conductor of 
changing length. 

To ensure that our system is behaving as expected, 
the theoretical limiting capacitance of an ideal cylindrical 
conductor must be compared to the measured capacitance of 
the cylindrical probe wires. The capacitance of an ideal 
cylindrical conductor is given by [Ref. 16] 




Eqn . 2.17 

where 
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e is the permittivity, e = Ks 0 ; 

K is the relative dielectric constant of a material; 

£ 0 is the permittivity of free space; 

1 is the length of the conductor, proportional to the 
height of the water; 

In is the natural logarithmic function; 
b is the radius of the outer charged surface; and 
a is the radius of the inner charged surface. 

The permittivity of free space is £ 0 = 8.854 x 10 - 12 

farad/meter . 



surprisingly low. A few inches of 0.2 mm diameter wire has a 
resistance to infinity of only 20 ohms. The capacitive 
impedance of our probe is many times higher, so that the 
water can be taken to be a perfect conductor. Electrical 
contact to the water is made in a two-wire measurement 
configuration through a second submerged wire at "infinity", 
many probe wire diameters away. 



unknown, unsteady, artifactual complex impedance between the 
probe and the bath at the very small length scales of the 



The resistivity of the water in our tank is 



A measurement of the capacitance gives the height, 

5c 

with the theoretical sensitivity ~ for the final capacitive 

ol 



test probes given by - 56.1 pF/m. 




The two-wire capacitive scheme suffers from an 
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order of a probe diameter. For example, the water meniscus 
slips and sticks to the probe because of a surface film that 
collects on the probes, instead of sliding freely as the 
water height oscillates. This film can itself add an unknown 
series resistance. And it is also possible that the 
conductivity of the water in the immediate vicinity of the 
probes is changing because of poorly understood 
electrochemistry at the metal-water-ion interface. For these 
reasons, we developed a four wire probe, based on the four 
wire impedance measuring technique commonly used in 
electronics and solid state physics. 

2 . Measurement Technique Models 

The physical electrical connections between an 
analyzer and a device under test (DUT) , water in our case, 
using a two-wire measurement technique can be modeled as 
unknown variable resistors in an equivalent electrical 
circuit as shown in Figure (2.1) . A constant current i is 
measured and applied to the DUT and the voltage v is 
measured across the series combination of the unknown lead 
impedance and the DUT impedance. 

I 1 

© 

& — CditQ — H 

Figure 2.1. Two wire measurement equivalent circuit. 
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Four wire, or four terminal pair configuration, 
measurements provide significant measuring advantage over a 
two wire measurement, for some situations. The HP-4192A and 
HP-4194A Impedance Analyzers are designed to make four wire 
measurements. The combination of equipment and measurement 
technique provides an accurate and easy measurement. 

A four wire measurement uses four connectors: High 
Current (H cur ) , High Potential (Hp 0t ) , Low Potential (Lp 0 t), 
and Low Current (L cur ) [Ref. 17]. The unknown variable 
impedances of the potential leads, as shown in the 
equivalent electrical circuit in Figure (2.2), are 
unimportant because there is no current in that portion of 
the circuit due to the high, nearly infinite, impedance of 
the voltmeter. Likewise, the unknown impedance of the 
current leads are unimportant because the current flowing 
through them is known, and all this current flows through 



the DUT. 



H 

cur 




L 

cur 



Figure 2.2. Four wire measurement equivalent circuit. 
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3. Four-Wire Conductance Technique 

Conductance is inversely proportional to the 
magnitude of the complex impedance. Conductance is 
proportional to changes in the water height. The conductance 
is obtained using the four-wire measurement technique by 
inverting the voltage difference in the water when measuring 
the complex impedance. 

The complex impedance Z of a device under test DUT 
is given by 



_ Pi _ 

Z " h w eff - i 



Eqn . 2.16 



where 



p is the complex resistivity of the conductor (water in 
our case) ; 

1 is the length of the conductor, proportional to the 
separation between the potential leads in the water; 

h is the height of the conductor on the wires; 

Weff I s the effective width of the conduction path 
between the wires; 

AV = Hp 0t - Lp 0t ; 

i is the current applied to the circuit. 

Figure (2.3) illustrates a four-wire resistance measurement 
in water. 
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Z y 



Water air 
interface 




Figure 2.3. Four-wire resistance measurement. 

Because it is based on the four-wire technique, our 
probe is insensitive to possible conductance changes of the 
water or to other contact impedances between the wire and 
the water on the "micro scale" illustrated in Figure (2.4). 
However, the four wire technique would be sensitive to 
changes in conductance on the "macro scale" between the 
probe wires. We do not know if this is a problem, or not. 



Macroscale 




Figure 2.4. Macro/micro scales about a four-wire sensor. 
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D. LEAST SQUARES FIT AND DEVIATION ANALYSIS TECHNIQUES 



The 


method 


of least 


squares for 


fitting 


data to a 


straight 


line 


can be 


accomplished 


on most 


hand-held 



scientific calculators today. The calculator automatically 
accumulates the necessary sums to provide slope, intercept, 
and linear correlation coefficient. Unfortunately, the 
correlation coefficient is not directly related to the 
uncertainty in the least squares values of the slope and 
intercept, which are the quantities of interest. [For a full 
development of the method of least squares fit, see 
reference 18] . 

The correlation coefficient r indicates whether there is 

causal relationship between the variables x and y, and 

therefore whether a linear fit to the data is justified. The 

correlation coefficient is defined by Taylor [Ref. 19] in 
terms of the covariance, a , and standard deviations a x and 

Gy . A more intuitive picture of the correlation can be 

developed by determining the slope m of a line of the data 

set when it is plotted y vs . x, and then when it is plotted 

x vs . y (which we will call m‘) . If there were no errors in 

the measurements and the measurement x and y obey a linear 

relation, then these two slopes would be a negative inverse 

of one another (i.e. m = -l/m‘) and r^ = 1, where r is 

defined as r = ; — If there is noise m the data, then r z < 



The quantitative significance of r is found when r is 

expressed in terms of probability. For any defined r Q , 
P N (|r|>r Q ) is the probability that N measurements of two 

uncorrelated variables would give a coefficient r as large 

as r Q . The probability that N measurements were derived from 

a random process involves the integration of a two 

dimensional Gaussian distribution [Ref. 20]. Generally it is 

more convenient to use a table which contains a value for r 

which would have to be exceeded for there to be a 
probability less than P N (|r|>r Q ) that the N measurements were 

generated by a random sample of observations taken from 
uncorrelated parent population. Two such tables are included 
in Taylor [Ref. 21] . 

For physical measurements, it is the error in the least 

square values of the slope and intercept which are of 

interest. Taylor gives expressions for the standard 
deviation in y, a , error in the slope, a m , and intercept, 

a^, in terms of the accumulated sums which are available in 

memory of most scientific calculators. I have reproduced 
them below: 



Oy 




Eqn. 2.16 
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n=l 



Z (5y n )2 = 



N 



n 



2 




Eqn . 2.17 




Eqn . 2.18 



a, 



m 



a. 



y 




Eqn. 2.19 



where n is a summation dummy index. 

Data sets were analyzed and errors in the fit parameters 
were determined using this technique throughout the thesis. 
Program code to determine the standard deviations in 
intercept, slope, and y using a HP-15C calculator is 
documented in Appendix A. The data was also analyzed using a 
Microsoft™ Excel file using equations (2.16) through 
(2.19) . 

It is also important to determine if there are random or 
systematic errors when examining experimental measurements. 
Random errors are revealed by repeating the measurements. 
Systematic errors cannot be revealed this way. Deviations in 
the least-square fit data are computed by subtracting the 
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regression calculated dependent value y c , from the measured 
dependent value y e [Ref. 22], 



Deviation = y e - y c . 

Eqn . 2.20 

Plotting the deviations helps reveal what type of errors are 
showing up in the measurement data . 

E. HEWLETT-PACKARD ANALYZER CROSS CORRELATION MEASUREMENTS 

"Correlation is a measure of similarity between two 
quantities" [Ref. 23]. This measurement technique allows you 
to measure a very small signal in the presence of noise. The 
dynamic signal analyzer can make this type of measurement 
well for very low frequency or single shot measurements, and 
provides a signal output in the time domain. 

There are two types of correlations: the auto- 

correlation function and the cross-correlation function. The 
results of the different correlations can be the same; 
however, the different correlations are applied to different 
types of systems. The auto-correlation function is a single 
channel measurement and the cross-correlation function is a 
two-channel measurement. 

The cross-correlation function compares two different 
signals as a function of the time shift between them [Ref. 
24] , and is given by 
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lim 

R xy ^ T ' ~ T->oo 



T 




0 



Eqn . 2.21 



where x is one input, y is another input, oo is infinity, 
and T is the time shift between the inputs. 

When the same signal is present in both channels in the 
presence of noise, the noise can be greatly removed using 
this technique. The correlation analysis is made when the 
dynamic signal analyzer is making a power spectral density 
measurement . 

Power spectral density measurements are made by the 
dynamic signal analyzer for each channel. A cross power 
spectral density is then computed internally, which is the 
product of the Fourier Transform of the two signals [Ref. 
25] , and is given by 



where / is the natural frequency and * indicates the complex 
conjugate of the function. 

The HP cross-correlation is then computed indirectly 
using the cross power spectrum [Ref. 24], and is defined as 



G xy<^ = S x (/)S y (/) 



Eqn. 2.22 



R,CyH> * F- 1 [G xy (/)l . 



Eqn. 2.23 
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F. GRAVITY WAVE PHASE SPEED AND WAVELENGTH MEASUREMENTS 

The phase speed of a sinusoidally driven gravity wave 
can be measured in a two-point measurement system using the 
output of the signal analyzer cross-correlation function and 
the distance between the sensor probes. This technique was 
used as a check of the capacitance probes. 

The phase speed of a wave is given by 



c = fk, 

Eqn 2.24 

and the theory of deep water gravity waves gives this speed 
as 



c 



2 ; if 



Eqn . 2.3 

Experimentally, the phase speed of a wave in the x-direction 
can be determined from the period of time T that a 
particular phase of the wave, such as a crest or trough, 
takes to move a given distance d. The phase speed is given 
by 



d 

c - T 

Eqn. 2.25 

where the distance d is the distance between the probes. 
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The cross-correlation output of a sinusoidal signal will 
also be a sinusoid, in the time domain. The period of a full 
or half cycle of a driven surface gravity wave can be 
determined from the cross-correlation measurement output. 
The period is determined using the method of least squares 
fit to the measured cross-correlation minimum or maximum 
data points . 

The phase speed is then a simple computation using 

equation (2.25). The driven surface gravity wave wavelength 

is determined using the measured period and computed phase 

speed. The natural frequency is the inverse of period, i.e. 
1 

/= ~. Substitute this into equation (2.24) and solve for 
wavelength to get 



X = cT. 

Eqn. 2.26 

G. FREE DECAY QUALITY FACTOR, Q MEASUREMENT 

The amplitude of a damped oscillator decreases 
exponentially with time. The amplitude may be represented by 

A ( t ) = A 0 e'P t , 

Eqn . 2.27 

where A 0 is the initial amplitude and P is called the 
temporal absorption coefficient. The length of time for a 
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damped oscillator to decay to 1/e of the initial amplitude, 
is called the decay time x, defined by 

1 2m 

T = P = V 

Eqn . 2.28 

where m is the mass and R m is the mechanical resistance 
[Ref. 26] . 

Substitute equation (2.28) into equation (2.27) and take 
the natural logarithm of the equation to get 

In [ A ( t ) ] = In [A ] - ; t. 

o x 

Eqn . 2.29 
1 

Using a change of variables, letting y = In [A(t)], m = 
and b = In [A 0 ] we get the equation of a line, 

y = mt + b. 

Eqn. 2.30 

For a continuous frequency system, such as the spectrum 
of a wind wave tank, the amplitude vs. time can be measured. 
Plotting In [A(t)] verses time, b and m can be determined 
using the method of least squares fit. The decay time can 
then be computed for any frequency f Q by 

1 

x = 

m 

Eqn . 2.31 
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Another method of reporting this decay is in terms of 
the quality factor, or "Q of the system" . The quality factor 
of a damped oscillator is defined [Ref. 27], 



Q = 



co 2 




Eqn . 2.32 

where ©2 and ©i ate average half power frequencies, above 
and below resonance © Q , respectively. 

I will not discuss the physics of a damped oscillator. 
This is discussed in "Fundamentals of Acoustics" [Ref. 28]. 
The Q can be expressed in terms of the mechanical constants 
of the system, given by 



Q 



V 
= Rm' 



Substituting the decay time from equation (2.28) 



Eqn 2.33 
gives 



Q = |© 0 t, or 
Q = nf Q x. 



Note that Q is a dimensionless quantity. 



Eqn. 2.34 
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III. EXPERIMENT 



A. INTRODUCTION 

The goal of our experimental program, of which this 
thesis is the first part, is to measure second-order 
(nonlinear) effects on deep water surface gravity waves. A 
very sensitive and linear measurement system is required so 
that transduction artifacts, such as non linearity, 
hysteresis, noise, etc., will be minimized. Interpretation 
of second-order effects will be much more reliable if we can 
ensure that we are using a linear, low noise, high 
sensitivity, wave height transducer system. 

There were two wave height transducer systems 
investigated in this thesis. The two-wire capacitance system 
was developed and tested until it was determined that the 
system was unsuitable for these applications. A four-wire 
conductance system was then developed and tested. The 
transducer systems development and testing are presented in 
parallel . 

The experimental apparatus can be divided into four 
major components, as shown in Figure (3.1). The first 
component is the wind-wave generating tank system. The tank 
system was originally built in the 1960 's by Professor 
Herman Medwin for use in experiments to measure sound 
scattering produced by wind driven waves on rough surfaces 
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[Ref. 29]. It was upgraded for these experiments. The second 
component is the sensor. Two sensor types were tested in 
these experiments. The third component is the associated 
electronics designed and built for each sensor. The final 
component is the off the shelf, general laboratory 
instrumentation . 




Figure 3.1. Experimental block diagram. 

This chapter will describe the general calibration 
techniques used, some of the design criteria and detailed 
construction of the experimental apparatus components (i.e., 
sensor, associated electronics, and the wind-wave tank 
system), and the specific calibration of all the apparatus. 
Finally, the techniques and experimental procedures used for 
the deep water surface gravity wave measurements will be 
detailed. 

The thesis will have several charts and tables that are 
referenced to the thesis lab books. The system will be lab 
book number, page number, represented by (B#P#). 

B. WIND-WAVE TANK SYSTEM 

1. Tank Facility 

This section will describe the capabilities and 
dimensions of the existing wind-wave tank and upgrades, the 
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first experimental component. The wind-wave tunnel section 
is 57 feet long, built primarily above ground, as shown in 
Figure (3.2). The tunnel is constructed of 3/4 inch plywood, 
fiber glassed for waterproofing, with 19 (2' x 4' x 0.5") (1 
x w x t) plate glass view ports. The water cross section of 
the tunnel is (4x3.5 feet) (w x d) . The wind cross section 
of the tunnel is (4x0.5 feet) (w x d) . The depths are, of 
course, only approximate, since they depend on the quantity 
of water placed in the tank. 

Top View 
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Figure 3.2. Wind-wave tank. 

The anechoic section is (10 x 10 x 9.5 feet) (1 x w x 
d) sunk into the ground, also shown in Figure (3.2). This 
section of the tank is constructed of six inch reinforced 
concrete and has (4" x 4" x 9.5') ( 1 x w x t) redwood wedges, 
used on the four sides and the bottom to absorb underwater 
sound in previous experiments. The top of this section has 
two removable 3/4 inch reinforced plywood covers for 
personnel safety. 
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2 . Small Wave Guide 



A paddle driven at large amplitude, can generate 
parametric waves that run transverse to the directly driven 
wave. These waves interfered with the dynamic calibration of 
the four wire probe. You can cutoff these unwanted waves 
using a wave guide. A small wave guide was constructed and 
placed in the main tank as shown in Figure (3.3). 
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Figure 3.3. Small wave guide. 

This wave guide has only two sides and a small 

paddle at one end. The top, bottom, and end opposite the 

paddle are open to the main tank. The minimum width of the 

wave guide was determined by the minimum width of the 

conductive probe, which was calibrated in this wave guide. 

The cutoff for this wave guide is 3.95 Hz. 

One side of this wave guide is a (2" x 6" x 10') ( t x 

w x 1) pine board. The other side is a (0.060" x 2’ x 8') (t 

x w x 1) Plexiglas plate. The Plexiglas is used so that the 

wave structure in the wave guide can be seen. The two sides 
are held together with six, eight inch long, ■g-16 bolts, 

3 

five, four inch Teflon spacers between the sides, and six, g - 
-16 hex nuts. 
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A (2 x 4 x 13 inch) (t x w x 1) pine board was used 

3 

for the paddle. An eight inch long, 'g-16 bolt went through a 
half inch hole that was drilled at one end of the paddle and 

3 

two (j x 5.5 x 16 inch) (t x w x 1) pine boards that were 

also bolted to the wave guide sides. A wave driver is 
connected to a 16 inch L-channel which was bolted to the top 
of the paddle. 

3. Wind Generation 

3 

The wind source consist of five pairs of ^ 

horsepower centrifugal fans on the top of a (39 x 48 x 41 
inch) (1 x w x h) housing. This housing is attached to the 
tunnel between locations one and two in Figure (3.2) . There 
are three power switches for the fan bank resulting in seven 
combinations of powered fans. 

The wind was measured approximately four feet from 
the wind-water interface, approximately in the center of the 
wind tunnel cross sectional area, at location two in Figure 
(3.2). Wind speed was measured using a hand held Turbo- 
meter™ manufactured by Davis Instruments of San Leandro, 
CA. The uncalibrated, measured wind speeds are reported in 
Table (3.1). 
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TABLE 3.1. WIND AND SURFACE DRAG SPEED AT LOCATION TWO 



(B2P55P131 ) . 





Wind 


Surface 


Fan 


Speed 


Drag Speed 


Combinations 


(m/s ± 0.1) 


(m/s ± 0.01) 


FI 


Too small to 


0.09 




measure 




F2 


2.7 


0.09 


F3 


2.9 


0.12 


FI + F2 


4.2 


0.12 


FI + F3 


6.5 


0.19 


F2 + F3 


6.0 


0.17 


FI + F2 + F3 


7.0 


0.22 



4 . Wave Driver 

The wave driver is a gravity wave source. Three wave 
drivers will be described in this section. 

a . Medwin Wave Driver 

The Medwin wave driver was developed and used in 
previous experiments in the tank system. The Medwin wave 
driver has a fixed amplitude and four drive frequencies. It 
consists of a forcing wedge and belt driven electro- 
mechanical plunging system. New wave drivers were developed 
for these experiments after initial tests in the tank 
determined that the Medwin wave driver was nonlinear itself. 
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The forcing wedge was a right triangle wedge 
that is pushed and pulled up and down in the water, Figure 

(3.4). The 1716 cubic inch right triangle wedge is 

3 3 

constructed using v inch plywood sides, braced with ^ inch 

plywood on both ends and in the center, and open on the top. 
There are one and a half inch guide strips on the outside of 
both ends of the wedge that are inserted into tracks built 
in the end of the tank just behind the fan housing. The 
tracks help direct the wedge movement up and down into the 
water. The entire wedge is treated with waterproof paint for 
submersion into water. Sand is uniformly placed inside the 
wedge to increase the weight in an attempt to increase the 
wedge stability. 

The plunging system consist of an electric 

motor, a fixed mechanical speed reducer, fly wheel, and two 

plunging rod assemblies shown in Figure (3.5). The electric 

motor is a Dayton Electric MFG. CO., model number 5K452C, 

half horsepower, 1725 RPM, 60 Hz, single phase, capacitive 

start motor. The motor is mounted on standard (1 x 1.5 x 10 

inch) motor mounts. The motor mounts are attached to a (11 x 
3 . 

8 x ~q inch) (1 x w x t) aluminum plate using four, one inch 

3 3 

long, "g-lO bolts. This plate is attached to a (13 x 22 x g - 

inch) (1 x w x t) aluminum bottom plate using an eight inch 

piano hinge and a four and a half inch long, ^-13 bolt. This 

construction allowed the position of the motor to be easily 

changed for belt adjustment and/or replacement. The large 
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bottom plate is shock mount to two (3.5 x 3.5 x 24.5 inch) (t 




x w x 1) redwood boards. These boards are clamped to 
aluminum L-beams that span the top of the tank, just behind 
the fan housing at location one in Figure (3.2). The 
plunging system is mounted above the water directly above 
the wedge . 



Figure 3.4. Medwin wave plunger. 

The mechanical speed reducer is an Ailing-Lander 

Co., Inc., model number CU-20, 20 to 1 ratio reducer. An 

aluminum ten inch diameter flywheel is attached to the ^ 

7 . 

inch diameter output drive axle. An eight inch long, g - inch 



diameter stainless steel rod, with two and a half inches of 
-13 threads is connected to the output flywheel. 
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Figure 3.5. Medwin drive system. 

1 

A Craftsman™, catalog number 9-1631, (j x 31 

inch) engine belt is used to connect the electric motor to 

the mechanical speed reducer. The speed reducer is mounted 

to the bottom plate also. The belt is used on two, four 

position, flywheels. Positioning the belt in any of four 

flywheel combinations determined which of the four output 

frequencies the system generated. 

The plunger rods assembly consisted of two (2 x 

2.5 x 1 inch) aluminum bearing housing, which are attached 
7 . 

to the o' inch diameter rod on the flywheel. The housings 
3 

have a •g-16 threaded hole that a 14.5 inch threaded rod may 

be screwed into and held with a lock nut at the housing. A 
seven inch turn buckle is attached to the other end of the 
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threaded rod. A two and a half inch eye bolt is attached to 

3 

the other end of the turn buckle. A six inch long ‘g-20 bolt 

with nuts goes through the eye bolt and the center pieces of 
plywood in the wedge, completing the plunger system. 

The non linearity was not quantified, but was 
readily apparent when watching the erratic plunging action 
and resulting waves of this driver. The wedge would not move 
up and down in a uniform motion. It erratically got caught 
in the guides due to the fact that the force from the 
plunger rods was not always normal to the wedge itself. This 
caused the wedge to be pushed or pulled from side to side as 
well as up and down. 

b. Yarber Wave Driver 

The Yarber wave driver is continuously variable 

in both frequency and amplitude. It has a forcing paddle and 

a direct drive electro-mechanical plunging system. 

3 

The paddle is a (12 x 44 x ^ inch) (w x 1 x t) 

3 

pine board with a (15 x 4 ^ inch) tab, which is just right 

of the center line as shown in Figure (3.6). The paddle is 

attached to two (2" x 4" x 4')(t xw x 1) boards with two, 

three and a half inch hinges. The two by fours are clamped 

to the end of the wave tank, using two six inch C-clamps. 

The paddle system consist of an electric motor, 

a direct drive speed reducer, and a single plunging rod 

assembly shown in Figure (3.7) . The motor is a General 

3 

Electric Co., model number 5KC45PG1140H, ^ horsepower, 1725 
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RPM, 60 Hz, single phase, capacitive start motor. The motor 
is connected to a Zero-Max© direct drive, model number QXl , 
CCW output, zero to 400 RPM continuously variable speed 
mechanical reducer. Full specifications are found in 
Reference 3 0 . 




Figure 3.6. Yarber wave plunger. 

The ten inch flywheel and plunger rod 

assembly were removed from Medwin's wave driver and 

connected to this system. The plunger rods are attached to 

3 

the flywheel using a new (12 x 2 ^ x 1 inch) (1 x w x t) 

slotted aluminum adapter. The slotted adapter provides a 
continuously variable adjustment of the driving amplitude. 
The plunging system was mounted to a one cubic foot wooden 
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box, that was clamped to a cement laboratory table. The 
cement table was placed as close to the tank as possible. 




Figure 3.7. Yarber drive system. 

The Yarber wave driver would change frequency at 
a rate of approximately 0.002 Hz per minute (B2P61) when the 
system ran continuously for over two minutes due to 
overheating in the mechanical speed adjuster. The amplitude 
and frequency were not adjustable during a measurement. 
This combination of things caused us to develop a third 
drive system. 

c. Keolian Wave Driver 

The Keolian wave driver is also continuously 
variable in both frequency and amplitude. It uses the same 
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forcing paddle assembly as the Yarber wave driver and an 
electric plunging system. 

The plunging system consist of a function 
generator, power amplifier, electric shaker, and a single 
plunging rod assembly shown in Figure (3.8). A HP-3314A 
function generator is connected to an Acoustic Power 
Systems, Inc. APS-114, dual-mode power amplifier. The 
amplifier is operated in the current mode. The amplifier is 
connected to an APS-113 Electro-SELS® shaker. A Kikusui, 
COS-5060A, 100 MHz oscilloscope is used to monitor the 

amplifier input and output. 




Figure 3.8. Keolian drive system. 

The plunging rod assembly consist of a 14.5 
3 

inch, g'-16 threaded rod with a one inch ball joint on each 

end. One ball joint was connected to a three inch aluminum 

1 

L-bracket, which was bolted to the shaker using inch long, 
■^■-20 bolts. The other ball joint was connected to a L- 

bracket on the paddle. 
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The wave amplitude and frequency are easily 
adjusted at any time using the synthesizer/function 
generator. Complex wave forms and noise can also be applied 
to the paddle using the amplifier and shaker. This will be 
important in latter experiments. 

5. Surface Drag Measurements 

Surface drag measurements were made at various wind 
speeds, approximately four feet from the wind insertion 
point, at location two in Figure (3.2). 

A special apparatus was constructed to hold the 
wires and associated electronics of the capacitive probe 
setup in the tank. We call this apparatus the Harp. The Harp 
is constructed of two 47.5 inch long aluminum U-channels, 
connected by two 30 inch long, 1/2 inch diameter aluminum 
rods. The upper U-channel is 2.2 5 inches wide by one inch 
high. The lower U-channel is 0.9 inches wide by 0.6 inches 
high. The wires were mounted using number one rubber bands 
attached to one end of the probe and the lower piece of the 
Harp. The upper ends were attached to a guitar adjuster, 
which was mounted to a 1/4 inch Plexiglas strip that was 
connected to the Harp. An electronics box was screwed onto 
the Harp. 

The Harp was put into the water with a single 
capacitance probe and a clear plastic ruler attached to the 
Harp with a three inch C-clamp. The' ruler and probe were 
separated by ten centimeters at the air-water interface and 
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the line formed was normal to the wind direction. Lint was 
placed on the water upwind of the ruler. The dust travel 
time between the two points was measured using a Casio Alarm 
Chrono G-Shock digital watch. The measured surface drag 
speeds are reported in Table (3.1). 

C. CALIBRATION TECHNIQUES 

1. Introduction 

There are two types of calibration techniques which 
were performed in this thesis. The apparatus and procedures 
for two calibrations will be described in this section. 

The first type of calibration was a static 
calibration which was used primarily during the probe 
selection. Static calibrations were also conducted to 
establish system baseline sensitivities. 

The second type of calibration was a dynamic 
calibration. The dynamic calibrations were the most 
important for the thesis because they set the sensitivity 
limits for the sensor system for all further measurements. 

2. Static Calibration 

The static calibration technique had three component 
parts as shown in Figure (3.9). The specific parts used as 
the water or electronics components could be changed, 
depending on the specific test conducted. A Starrett Co., 
No. 254, 26 inch master vernier height gage was used in all 
static calibrations. 
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The static calibration technique is simple. The 
sensor or sensor system was somehow attached to the height 
gage. The necessary measurement electronics were connected 
to the sensor or sensor system. The sensors were immersed in 
water. The sensor parameters or sensor system electrical 
outputs were measured, verses the measured relative 
immersion determined using the height gage. A sensitivity 
was then computed, in the form 

A output 
A height ' 

Eqn 3 . 1 




Figure 3.9. Static calibration components. 

3. Dynamic Calibration 

The dynamic calibration techniques also have three 
component parts as shown in Figure (3.10). All of these 
components could be changed depending on the specific test 
conducted, and components could be combined during certain 
tests. The dynamic calibrations were limited to sensor 
system tests. 

There were two movement possibilities for the 
dynamic calibrations. The sensor system could be moved in 

42 




the water or the water could be moved past the sensor 



system. Dynamic calibrations for both movements were 
attempted . 




Figure 3.10. Dynamic calibration components. 

A subjective decision was made that the best dynamic 
calibration was an dynamic calibration for a constant 
amplitude traveling water wave, moving past the sensor 
system. This decision was made because the traveling wave is 
a gravity wave. Gravity waves and wind-water waves were what 
the sensor system was designed to measure. The other dynamic 
calibration techniques attempted are described fully in 
Appendix B. 

The sensor system was placed in the small wave guide 
describe in chapter III-B-2-b. A static calibration was 
conducted to establish the baseline sensitivity. A 12 inch 
long, 0.031 inch diameter brass rod was scotch taped to a 
Starrett Co., No. A LH, three axis micrometer. The 
micrometer was bolted to the main tank approximately three 
feet from the paddle and directly above the small wave 
guide . 
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The brass rod was raised 0.354 cm above the 
unperturbed water. This would provide a visual point with a 
0.25 cm RMS maximum for the traveling waves. 

The Keolian wave driver was attached to the paddle 
of the small wave guide. The drive amplitude was adjusted so 
that the traveling wave amplitude was constant, ± 0.1 cm, as 
determined by the peaks of the wave just touching the bottom 
of the brass rod for several frequencies below 4.0 Hz. The 
small wave guide cutoff for the paddle generated parametric 
waves was 3.95 Hz which limited the frequency of the 
calibration. The sensor system electrical output for a 
constant amplitude traveling wave was measured, verses 
applied drive frequency. A sensitivity was then computed, in 
the form 

A output 
A frequency’ 

Eqn 3 . 2 

D. WAVE HEIGHT PROBES 

The second experimental component is the wave height 
sensor probe that we need to make the surface gravity wave 
measurements. The goal is to develop a height sensor that is 
linear, very sensitive, and noise immune. A variety of wave 
probes have been used by experimenters using changes in 
capacitance and resistance [Ref. 5, 7] . Two sensor probes 
were tested in this investigation. The constraints in the 
probe selection, test procedures, results, and final probe 
selection will be discussed in this section. 
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1. Capacitance Probe 



The capacitance probe investigated in this thesis 
utilized the two-wire measurement technique described 
earlier. The variable capacitance measured is determined by 
the probe immersion into water. The capacitance probes were 
all labeled and marked T-(#) so that repeated measurements 
could be made for analysis for each probe. The probes were 
placed on the Harp for all capacitance test conducted in the 
wind-wave tank. 

There were two major categories of constraints on 
the capacitance probe selection. Optimizing the various 
electrical and hydrodynamic constraints was an interesting 
challenge . 

There are three electrical constraints on the probe 
when it is in the water. The first constraint is that the 
resistance between the probe and water should be nearly 
infinite . 

The second constraint is that the capacitance per 
unit length should be as high as possible to maximize the 
sensor sensitivity. From equation (2.15), the capacitance 
increases as the radius of the outer charged surface 
approaches the radius of the inner charged surface. We 
wanted as thin a dielectric as possible. Meeting this 
constraint would establish an "ideal" electrical capacitive 
component . 
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The third constraint is that the capacitance 
sensitivity, should be constant over the length of the 
cylinder, to remove unwanted systematic errors. This 
required a uniform dielectric. 

There were also two interrelated hydrodynamic 
constraints on the probe that also require optimization. The 
first constraint was that the probe drag needed to be 
minimized to reduce the induced turbulence in the flow 
field, caused by the probes. This also required a thin 
conductor and dielectric. Second, the effect of probe 
wetting needed to be minimized so that the change of 
capacitance closely followed the water height, reducing 
hysterisis . 

The optimization of the hydrodynamic constraints had 
to be accomplished on two time scales. The period of a 
forced surface gravity wave is the first time frame. This 
time frame is important, as it sets the low frequency limit 
of the sensor. The period of an experimental data run is the 
second time frame. The capacitance sensitivity, must remain 
constant during this time frame to prevent systematic errors 
in the measurements. 

Finding the right sensor seemed like it would be an 
easy job. First we needed to find an electrical conductor 
that had nearly infinite resistance, that could be used as a 
conductor/dielectric/water capacitor. 
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Determining what wire should be used was the first 
order of business. The first test was to find a linear and 
highly sensitive probe. The static calibration technique 
discussed in chapter III-C-2 was used to measure the probe 
sensitivity. The sensor probes were submerged in water taken 
from the tank and contained in a (2.25 x 17.5 inch) (diameter 
x h) , 1000 ml, Pyrex™ graduated cylinder. The capacitance 
of a submerged probe was measured using a HP 4192A LF 
Impedance Analyzer. 

Measurements were made with the HP-4192A using a 
30.0 kHz, 1.0 volt drive signal. The 30.0 kHz signal is used 
because comparable frequencies will be used in latter 
experiments. A 5.0 ounce plum bob is used to straighten the 
probe and insure that it is perpendicular to the air-water 
interface. Figure (3.11) is a diagram of the experimental 
apparatus . 

HP-4192A 



Starrett Co . 
Master Vernier 
Height Gage 
No. 254 



■t*. 




Sensor 

Ground 





Pyrex 1000 ml 
Graduated Cy 1 inder 

Plumb Bob 
-Tank Water 



Figure 3.11. Static capacitance sensitivity experimental 

setup . 
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The first capacitance probes that were tested were 
all made from magnet wire. Several coatings were tested on 
the magnet wire. There were eight probes built, using 22 
gauge magnet wire manufactured by Belden and Consolidated. 
The ends of half the probes were dipped in blue Testors 
model paint. The ends of the other half of the probes were 
coated in BIPAX, TRA-BOND, BA-2106, 5-minute epoxy. Half of 
the probes were coated in Bulls Eye Shellac, while the 
remaining probes only had their manufacturer's insulating 
coating. Neither the capacitance nor the dissipation were 
constant, at a constant immersion depth, for a length of 
time comparable to an expected gravity wave measurement. 

A second capacitance probe was built, but not fully 

tested because the measurement technique was changed. A 33 

1 

inch long, 0.160 inch ID., and ^ inch OD . glass capillary 

tube formed the dielectric for this probe. A Malco S-50 
Micro dot connector, P.N. 031-0001-0001, was soldered to 22 
gauge tin copper wire (Aztec Wire & Cable) . One end of the 
glass was melted closed and the wire was place in the 
capillary tube. The connector end of the wire was attached 
to the glass with several layers of shrink tubing. 

The final wire tested and selected for use was Omega 
Engineering, Inc., TF-010, 10 mil diameter, Teflon coated 

copper thermocouple wire. The Teflon has a dielectric 
constant K = 2.0 to 2.1 [Ref. 31] . It was selected for the 
sensor probe because its capacitance was constant over time. 
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Figure (3.12) shows measurement samples of capacitance 

verses time for shellacked magnet wire (T-2) and Omega 

thermocouple wire (T-15) . The shellacked magnet wire 

pF 

capacitance changed at a rate of 4.5 The Omega 

pF 

thermocouple wire capacitance changed at a rate of 0 . 013jhj“. 




Figure 3.12. Capacitance of immersed probe wire vs. Time 

(B1P36) . 

The immersed end of the Teflon wire was looped and 

tied to the plum bob or other tensioning device. The knot 

was covered with Polymeric Systems, Inc., PSI 601 Clear, 

1 

100% Silicone, RTV Adhesive Sealant. Approximately •j inch of 

electrical heat shrink tubing was shrunk over the knot and 
RTV. The "dry” end of the probe was tinned and soldered into 
a male, 0.062 inch, terminal (Digi-Key part number WM1000). 
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Test were performed in tank water and in mercury 
using the static equipment setup shown in Figure (3.4). The 
mercury was held in a different container for safety 
reasons. The mercury test should give results as close to 
theory as experimentally possible. 

The theoretical capacitance per unit length for this 
geometry using equation (2.15) is Cj-heory = 2.40 ± 0.01 
pF/cm (B2P73). The measured capacitance sensitivity was 2.2 
± 0.1 pF/cm in mercury and 2.4 ± 0.1 pF/cm in tank water, as 
shown in Figure (3.13), and was linear as a function of 
immersion in both cases. 
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Figure 3.13 Capacitance of Teflon wire in different 
media vs. Relative immersion (B2P22/102). 

A wire capacitance static calibration was performed in 
the wind-wave tanks using precision blocks for the height 
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gage. The blocks were 1.7 ± 0.1 cm high. The measured 
capacitance sensitivity for this test was 2.2 ± 0.1 pF/cm, 
as shown in Figure (3.14), and was linear as a function of 
immersion . 

A final static capacitance sensitivity test was 
performed in the large wind-wave tank using a HP4194-A 
Impedance/Gain Phase Analyzer. The four wire measurement 
technique was used to measure resistance R and reactance X. 
The capacitance probe was on one side of the DUT, tank 
water, and the other side was the Harp or a paper clip. 
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Figure 3.14. Capacitance vs. Relative immersion for T-15 
using precision blocks (B2P108) . 

The capacitance can be determined for the measured 
reactance, as a function of frequency, [Ref. 31] by 



C (/) 



1 

271 /X' 



Eqn 3 . 3 
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The probe wire was mounted on the Harp and immersed 
approximately 25.0 ± 0.5 cm. The theoretical capacitance for 
a cylinder of this length is 60.0 pF . The measured 
capacitance was 62.5 ± 0.1 pF at 30 kHz using the paper 

clip, and 66.5 ± 0.1 pF using the Harp, as shown in Figures 
(3.15 and 3.16). 

The static capacitance sensitivities for the various 
measurement techniques are shown in Table (3.2). The 
measured sensitivities are in general good agreement with 
theory . 
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Figure 3.15. R & X vs. / for Teflon wire (T-15) in the wind- 
wave tank using a paper clip (B2P93) . 
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Figure 3.16. R & X vs . / for Teflon wire (T-15) in the wind- 
wave tank using the Harp (B2P93). 
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TABLE 3.2. STATIC CAPACITANCE SENSITIVITIES FOR ALL 



MEASUREMENT TECHNIQUES . 



Method 


Sensitivity 
(±0.1 pF/cm) 


Cylindrical Theory for TF-101 


2.4 


T-13 immersed in graduated 
cylinder with tank water at 30 
kHz (B2P23 ) 


2.2 


T-15 immersed in graduated 
cylinder with tank water at 30 
kHz (B2P23 ) 


2.3 


Omega TF-101 in Mercury 

(B2P102) 


2.2 


T-15 immersed on Harp in wind- 
wave tank using precision 
blocks (B2P108) 


2.2 


T-15 immersed on Harp in wind- 
wave tank at 3 0 kHz, Harp 

qround (B2P93) 


2.7 


T-15 immersed on Harp in wind- 
wave tank at 30 kHz, paper clip 
qround (B2P93) 


2.4 



54 



2 . Conductance Probe 



The idea of a four wire measurement was motivated by 
continued electrical grounding problems in the tank when 
using the capacitance probe and the inability to accurately 
measure the electrical properties of water. The four wire 
technique was described in chapter II -D. It was not obvious 
however, what quantity should be measured. This section will 
describe how we decided to measured conductance and the 
selection process for the best conductance probe. 

The type of measurement was the first item to be 
determined. The static calibration technique discussed in 
chapter III-C-2 was used for these test. Tank water was 
placed in either a 1000 ml Pyrex™ graduated cylinder or a 
four liter desiccator. The BNC cables were connected to a 
HP-4194A Impedance/Gain-Phase Analyzer in the four wire 
measurement mode. A Think Jet printer was connected to the 
analyzer by HBIB bus for data recording. 

A (11.5 x 43.0 cm) (1 x w) Vector perforated printed 

circuit board had 20, number 43 holes, drilled 0.125 inches 

1 

from one end. The holes were 0.0125 inches apart so that 

inch, 2-56 stainless steel machine screws could be placed in 
the board. These machine screws became the restraining 
mechanism and electrical connectors for the various samples. 
Two BNC cables were tied to the Vector board. Two Kings, KC- 
79-35, E-Z mini hooks were then attached to the machine 
screws and BNC cables. The Vector board was clamped to a 
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Starret Co., No. 254, master vernier height gage with a one 
inch C-clamp. The experimental setup is shown is Figure 
(3.17). 



Think Jet Printer l ~~1 




Figure 3.17. Conductance experimental test setup. 

Resistance R, reactance X, conductance G, and 
susceptance B, were all measured as a function of wire 
immersion. A 30 kHz, one volt signal was used to drive the 
four 0.0195 inch diameter brass wires that were 0.125 inches 
apart. Brass was used only because it was most accessible 
at the time. Conductance proved to be linear as shown in 
Figure (3.18), which is the type of measurement we want. 

The next test, using the same equipment setup, was 
conducted to determine what kind of metal and sensor 
geometry was best, if there was a frequency dependence, and 
if there was any dependence on drive signal or time. 
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Using a swept frequency in two ranges, 100 Hz to 50 
kHz and 100 Hz to 2 MHz, R, X, G, and B were measured as a 
function of depth for each probe. This provided a check for 
any frequency dependence in the measurement. Typical test 
measurement data is shown in Figures (3.19 & 3.20). The x 
and y axis are both linear in all cases. There is no 
frequency dependence evident. 



Conductance (mS) or 
RaaMtanca (Ohm) 



♦ Conductance G 
O Reactance R 



Relative Unmartion<cm| 



Figure 3.18. Brass wire G & R vs . Relative immersion 

(B2P118) . 

Brass, gold, platinum, and stainless steel were 
tested in a V shaped or cylindrical geometry. The brass were 
(10.4 cm x 0.288" x 0.005") (1 x w x t) V shaped strips and 
0.195 inch diameter, 10.4 cm long wires. The gold was 0.20 
inch diameter, 8.9 cm long wires. The platinum was (16.5 x 
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0.68 x 0.032 cm) (1 x w x t) V shaped strips. The stainless 
steel was 0.20 inch diameter, 10.2 cm long wires. 
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Figure 3.19. R & X for a stainless steel wire at a 6.0 cm 
relative immersion vs. Frequency (100-50,000 Hz). 
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Figure 3.20. G & B for a stainless steel wire at a 6.0 cm 
relative immersion vs. Frequency (100-50,000 Hz). 
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The stainless steel wire had the highest sensitivity 

mS 

of the metals tested, at 8.4 — . The brass wire had a higher 

sensitivity than the brass strip. We decided to use 
stainless steel wire because of its high sensitivity. A plot 
of the measurement data is shown in Figure (3.21). 
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Figure 3.21. Conductance of four metals and two geometry's 
vs. Relative immersion (B2P114-119 ) . 

Figure (3.22) shows that there is no dependence on 

oscillator drive for brass or stainless steel wires, at a 

constant immersion depth of approximately seven centimeters. 

The final test for the conductance probe wire 

selection was to determine if there was a time dependence 

for the stainless steel wire. There is a very small 
dependence of 0.01~^~, shown in Figure (3.23). 



60 



140 _ 



100 -- 



80 -- 

Conductance ImS) 



40 -- 



20 __ 



♦ EV as* W»e 
■ Stainless Steel Wie 



-i 1 1 ! i 1 

0 0.2 0.4 0.6 0 8 1 

Ovcitaior <krv# (V) 



Figure 3.22. Conductance of two metals vs. Oscillator drive 

(B2P148) . 
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Figure 3.23. Conductance of stainless steel wire vs. Time 

(B2P149) . 
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A five wire configuration of stainless steel wire 
was tested and selected for use in the conductance method. 
This configuration was selected for the probe because two 
separate two wire measurements could be made simultaneously 
at approximately the same location. It was anticipated that 
the two outputs might be used in a gradient measurement. 

The probe has five, 0.020 inch diameter, five inch 
long stainless steel wires that are held at the top by a 
(1.5 x 1 inch) Vector board, and held at the bottom by a 
(1.5 x 3/8 inch) double sided copper printed circuit board 
(PCB). Five 0.020 inch diameter holes were drilled 1/4 inch 
apart, in a line, on both the PCB and the Vector board. 

A Dremel® Moto-tool®, model 395, was used to cut 
electrical pads for each hole on both sides of the PCB. The 
stainless steel wires were inserted in the holes of each 
board. The wires were then soldered on each side of the PCB. 
A 1/4 inch diameter semi-circular loop was formed at the top 
end of the wire, and one end was soldered into a Vector 
board hole along with a 20 inch long, 22 gage, wire used to 
make electrical connection to the associated electronics. 
The loop was used to adjust the tension of each individual 
stainless steel wire. 

A 2.5 inch diameter "J" was formed on one end of a 
20 inch long, 1/8 inch diameter, brass rod. The short side 
of the J was inserted into a 1/16 inch hole in an 
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electrically isolated pad on the PCB and soldered on both 
sides . 

A two inch right sided triangular sheet of brass was 
soldered to the other end of the brass rod. A two inch long, 
8-32, brass bolt was soldered to the opposite side of the 
triangular sheet. The screw was inserted into a 1/4 thick 
piece of Plexiglas which supported the Vector board. The 
Plexiglas gave the top of the probe some rigidity. A 8-32 
brass nut was used on the bolt to adjust the tension on all 
five wires. 

The brass rod was bolted to a 3/4 inch diameter PVC 
pipe. The 2 0 inch wires were fed through the PVC. The PVC 
could then be clamped and the conductance electronics could 
be mounted to the other end of the PVC. The five wire probe 
is shown in Figure (3.24) . 
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E. MEASUREMENT SYSTEM ELECTRONICS 



The third experimental component is the associated 
electronics, which will convert whatever the sensor is 
measuring into a usable electric signal. This component must 
also be linear, sensitive, and immune to noise. The 
constraints and design for the electronics, test procedures, 
results, and final electronics selection will be discussed 
in this section. 

1. Capacitance System Electronics 

The capacitance system electronics is a two channel, 
capacitance bridge, op-amp circuit. The purpose of this 
circuit is to provide an electrical output for the sensor 
probe by converting the capacitance changes of the probe 
wire into voltage changes which could be input to 
conventional electronic signal processing instrumentation. 
The two sections of the circuit must combine to give a very 
sensitive and linear output which represents the height of 
the wave on the sensor probe. 

The first section of the circuits are the 
transformer fed capacitance bridge, as shown in Figure 
(3.25), which provides a means to detect small signals. A 
lock-in amplifier oscillator output is fed into the bridge 
via transformer T1 . The transformers are manufactured by 
ADC, model number TF4RX17YY. They have a frequency range 
from 40 Hz to 100 kHz. 
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The sensor probe and water ground form one leg of 
the bridge, variable capacitor Cpj-o^g. The other leg has 
three capacitors in parallel, one fixed value Cl, and two 
variable value. There is a rough adjust C5 and a fine adjust 
C4 . The variable capacitors are used to electrically null 
the bridge when the sensor is in the tank water and the 
water is in an equilibrium state. The bridge can also be 
nulled when a dummy load capacitor is inserted into the 
circuit instead of the tank water. 

The typical equilibrium capacitance of the sensor 
probe is approximately 60 pF. The other leg of the 
capacitance bridge is the sum of the fixed, rough adjust, 
and fine adjust capacitors. The rough adjust capacitor has a 
value of 27 ± 13 pF. The fine adjust capacitor has a value 
of 10 ± 5 pF. The fixed capacitor in the adjustment leg is 
then 



c • = c - c - c 

fixed sensor rough fine* 



Eqn. 3.4 



The approximate value for Cf-j_ xec } is 23 pF. 

The second section of the circuit provides output 
impedance matching and signal gain using an op-amp. The OPA- 
111 op-amp is used with a gain adjust capacitor C6. The OPA- 
111 provides a typical open-loop gain of 120 dB and typical 
common-mode rejection of 100 dB. Complete characteristics 
are found in the manufacture's specification sheet. [Ref. 



32 ] 
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The two channels are identical except for the ground 
used after the transformer. The different grounds are used 
for circuit isolation and to minimize the channel cross- 
talk. The two channel outputs are matched and maximized 
using the gain adjust capacitor C6, in the op-amp feedback 
circuit . 

The electronics had one linearity test. Linearity is 
the most important feature required for the dynamic wave 
height measurement system. The principle of superposition 
can be used if the sensor probe and electronic system 
components each have linear responses. The electronic system 
is measured using "dry" techniques. 

The "dry" calibration of the electronic system is a 
lab bench experiment using mylar capacitors as a substitute 
for probe wire water capacitance. The values of these 
capacitors were all measured at room temperature using a HP- 
4192A Impedance Analyzer. The nine capacitors used for the 
dummy load test were measured and the results of this test 
are shown in Table (3.3). 

A EG&G Princeton Applied Research, model 5210, 
Lock-In Amplifier was used as a signal source and detector. 
A 31.3 kHz, 1.999 volt signal was applied to the input of 
the electronic system. When both channels are being used, 
two lock-in amplifiers were used. 
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Figure 3.25. Capacitance system electronics circuit. 
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TABLE 3.3. DUMMY LOAD CAPACITANCE VALUES (B2P16). 




Different prime numbers are used for the frequencies 
of both lock-in amplifiers. We used 28.3 and 31.3 kHz 
throughout this experiment [Ref. 33]. 

The output of the electronic system was applied to 
the input of the lock-in amplifier. The typical lock-in 
amplifier setup configuration is shown in Table (3.4). A 
Kikisui, COS 5060, dual channel oscilloscope was used to 
monitor the electronics system input and output signals. 
Figure (3.26) shows a block diagram of the experimental 
setup . 

TABLE 3.4. LOCK-IN AMPLIFIER SETUP CONFIGURATION. 



Control Button 


Setting 


Sensitivity 


3 00 mV 


Filters 


Flat & 
Tracking 


T 


30 msec 


Slope 


12 dB 


Dynamic 

Resolution 


Normal 


Display 


XYV 


Reference 


Internal 
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Figure 3.26. Dry electronic system experimental setup. 

The first step of the dry calibration is to zero the 
bridge. The variable capacitor C4 and C5 are adjusted to 
zero the output voltage when a 54.0 pF dummy load capacitor 
is placed in each circuit. 

The second step in the dry calibration is to match 
and maximize both channel outputs. Variable capacitor C6 is 
adjusted on both channels when a 10 pF dummy load capacitor 
is moved back and forth between the two channels. The output 
from a 450 pF range of capacitors is then measured to test 
the linearity of the electronic system. 

Channel one sensitivity was measured at 3.21 ± 0.02 

mV mV 

Channel two sensitivity was measured at 3.41 ± 0.03 ^r. 

The linearity for each channel is illustrated in Figure 
(3.27) . 
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Figure 3.27. Capacitance system electronics output voltage 
vs. Dummy load capacitance (B2P17). 

2. Conductance System Electronics 

The five wire probe illustrated in Figure (3.24) is 



connected to 


two lock-in amplifiers to make 


a 


four wire 


water surface 


height 


resistance measurement, 


as 


shown 


in 


Figure (3.28) . 












Conductance 


is inversely proportional 


to 


the 


magnitude of 


the 


complex impedance, or 


AV, 


but 


is 



proportional to the water height on the probe as described 
earlier in Chapter II (C) (3). The conductance system 
electronics is a two channel division circuit. The purpose 
of this circuit is to invert and balance the two lock-in 
amplifier output signals. The two channels are identical 
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because we thought we might perform a gradient measurement 
using the difference outputs. Gradient measurements were not 
completed . 




22345 Oscilloscope 

Probe 

Figure 3.28. Five wire probe electronics setup. 

The output voltage of a lock-in amplifier used in 
the differential mode, provides the input to an Analog 
Devices AD 534 Internally Trimmed Precision IC Multiplier, 
as shown in Figure (3.29). The output of the divider is 
determined by [Ref. 34]: 



10 Z2 

v out = + y l' 

Eqn 3 . 5 

where x^ is the signal input, y^ is the summing input, and 
Z 2 is the numerator input. 
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Resistors through R 5 form a voltage divider 
network which determines the overall gain of both channels. 
Resistor R 3 is adjusted to equalize the two channel outputs 
for the same input signal. The circuit is powered by a KEPCO 
Multiple Output Power Supply, MPS 62 OM. 
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Figure 3.29. Conductance system electronics circuit. 



72 




F. WIND-WAVE TANK EXPERIMENTS 



1. Introduction 

There were five measurements made in the wave-wind 
tank. They were the wave power spectrum, surface gravity 
wave phase speed, surface gravity wave wavelength, and wave- 
wind tank system free decay Q and impedance measurements. 
This section will detail the apparatus and techniques used 
for each measurement . 

2 . Power Spectrum Measurement 

All the wave-wind tank measurements except the tank 
impedance were based on the output of the power spectrum 
measurements as described in earlier sections. A signal was 
detected using either the capacitance or conductive 
measuring technique and then power spectrum measurements 
were made using a HP-3562A Dynamic Signal Analyzer. 
Measurements were made at various positions in the tank, and 
various wave-wind states. The basic states are: 

a) baseline noise, 

b) various wind speeds, 

c) various driven surface gravity wave speeds, and 

d) combinations of wind speeds and driven surface 
gravity wave speeds. 

Data was recorded using a HP-7470A plotter or stored 
on disk using a HP-9133 hard disk. To store data on disk, 
the HP-3562A front panel control sequence in Appendix C was 
used. A HP-54501A, 100 MHz digitizing oscilloscope was used 
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to monitor the different electrical output signals. The 



experimental electronics setup is shown in Figure (3.30). 
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Figure 3.30. Measurement electronics setup. 

Power spectrum measurements were taken with the HP- 
3562A using power spectrum measurement and display modes. 
Typical HP-3562A setup configurations are shown in Figure 
(3.31) . 
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Figure 3.31. Typical HP-3562A power spectrum setup 

configuration . 

3. Phase Speed and Wavelength Measurements 

The phase speed and wavelength were measured using 
the cross-correlation output of the HP-3562A Dynamic 
Analyzer. Two signals had to be detected to make the phase 
speed measurement as described in chapter II-F using the 
experimental setup shown in Figure (3.30). 

The HP-3562A computed and displayed the cross- 
correlation of the two signals. Typical HP-3562A setup 
configurations are shown in Figure (3.32). The Phase speed 
was then computed using the method of least squares fit on 
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the data output of the cross-correlation measurement. The 
distance between the two probes had to be measured, then the 
wavelength could be computed as described in chapter II-B. 
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Figure 3.32. Typical HP-3562A cross correlation setup 

configuration . 

4 . Free Decay Q Measurement 

The free decay Q was calculated using the power 
spectrum measurement output of the HP-3562A over a period of 
ten minutes. The ten spectra were measured one minute apart 
and stored on disk for each data set using the experimental 
setup shown in Figure (3.30), as described in chapter II-H. 



76 



The entire data set was latter plotted on the same 
paper. The data required for any measured frequency could be 
extracted from the data on the disk files, and the Q could 
be calculated. 

5 . Tank Electrical Impedance Measurement 

A self contained mini experiment was performed to 
answer the question of why the lock-in amplifier oscillator 
was being loaded down by 5-50% during various capacitance 
measurements. The tank impedance, via the water and 
plumbing, returning to the rack ground was tested. 

The lock-in amplifier output voltage was measured on 
an oscilloscope with and without the tank impedance in the 
circuit using the equipment setup shown in Figure (3.33) . 
The tank impedance is then calculated using the measured 
voltage at any given frequency, by 



Z tank 



v closed 

v 

open 



Z 



LI' 



Eqn . 3.6 

where Zj_j is the lock-in amplifier output impedance. The 
lock-in amplifier has a 600 ohm output impedance which was 
verified with a 100 ohm resistor. 
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COS-5060A 




Figure 3.33. Tank electrical impedance experimental setup. 

Voltage output for the closed circuit as a function 
of frequency is plotted and shown in Figure (3.34). The open 
circuit voltage was constant at 5.7 V pp . The measured tank 
impedance at 20 kHz is Z ( 2 0 kHz) = 21.5 ohms. 
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Figure 3.34. Lock-in amplifier output, voltage vs. Frequency 

for tank (B2P107). 
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IV. RESULTS 



A. STATIC CALIBRATIONS 

The static calibrations were certainly the easiest 
measurement of the experiment. They were easy because there 
was no complicated equipment, the technique described in 
chapter III-C-2 is very easy to follow, and we used the 
static calibration technique so many times. 

1. Capacitance Static Calibrations 

All of the capacitance probe static calibrations 
were in good agreement with the theoretical capacitance 
sensitivity of 2.4 ± 0.1 pF/cm as shown in Table (3.2). 

A theoretical capacitance system static voltage 
sensitivity could be determined by combining either the 
theoretical capacitance sensitivity or the measured Teflon 
wire capacitance sensitivity with the measured electrical 
sensitivities of each channel of the system electronics. The 
theoretical voltage sensitivities for each channel are shown 
in Table (4.1) using the results from the previous chapter. 

Capacitance system static voltage sensitivities were 
measured using the static calibration technique with water 
contained in the Pyrex graduated cylinder and the 
electronics shown in Figure (3.26). The lock-in amplifier 
oscillator was set at 30 kHz for these measurements. A plot 
of the measurement data is shown in Figure (4.1). 
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TABLE 4.1. THEORETICAL CAPACITANCE SYSTEM STATIC CALIBRATION 



ITEM 


SENSITIVITIES 


dC ( qf ) 
dH Vcny 

±0.1 


dV fmV"\ 
dC IpFj 


dV fmV\ 
dH vcny 

±0.1 


Theoretical Teflon 
TF -101 


2.4 






T-13 in Graduated 
Cylinder (GC) (B2P23) 


2.2 






T-15 in Graduated 
Cylinder (GC) (B2P23) 


2.3 






T-15 in tank (T) 
(B2P108) 


2.2 






Chi (B2P18) 




3 . 21±0 . 02 




Ch2 (B2P19 ) 




3 . 41±0 . 03 




TF-101 in Chi 






7.7 


T-13 (GC) in Chi 






7.1 


T-15 (GC) in Chi 






7.4 


T-15 (T) in Chi 






7.1 


TF-101 in Ch2 






8.2 


T-13 (GC) in Ch2 






7 . 5 


T-15 (GC) in Ch2 






<1 

00 


T-15 (T) in Ch2 






7.5 
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Figure 4.1. Capacitance system electrical output voltage vs. 

Relative immersion in the graduated cylinder (B2P24). 

Capacitance system static calibrations were also 
conducted using the Cone Head apparatus described in 
Appendix B. Water was added to the Cone Head before each 
measurement during this test, using Omega thermocouple wire 
probe T-13. A plot of the measurement data is shown in 
Figure (4.2) . 

The capacitance system static voltage sensitivities 
for channel one measured in the graduated cylinder are in 
good agreement with theory. It is unknown why the 
capacitance system voltage sensitivities for channel two 
measured in the graduated cylinder were not in better 
agreement with theory. 
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Figure 4.2. Capacitance system electrical output voltage vs. 

Relative immersion in the Cone Head apparatus (B3P2) . 

The system voltage sensitivities measured using the 
Cone Head apparatus are not in good agreement with theory. 
This is probably because of the apparatus and procedure 
itself. We could accurately measure the water that was being 
put into the Cone Head system. We could not accurately 
measure the amount of water that was actually getting into 
the changing water column. Water would stick on the sides of 
the graduated cylinder and the wire. 

The theoretical and measured sensitivities for all 
the test, and the percent differences are shown in Table 
(4.2). 
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TABLE 4.2. CAPACITANCE SYSTEM STATIC SENSITIVITIES (B2P29). 



System Combination 


Theoretical 

Sensitivity 
dV /in Vi „ 

dH vcmj ±0,1 


Measured 

Sensitivity 
dV /in Vi „ 

dH Vcmj ±0 * 1 


A% 


T-13 & Chi in 








graduated cylinder 
(GC) 


7.1 


7.3 


2.8 


T-13 & Chi in Cone 
Head (CH) 


7.1 


6.45 


-9.86 


T-15 & Chi (GC) 


7.4 


7.4 


0.0 


T-13 Sc Ch2 (GC) 


7.5 


8.1 


8.0 


T-13 Sc Ch2 (CH) 


7.5 


8.5 


13.33 


T-15 Sc Ch2 (GC) 


7.8 


8.1 


3.85 



2. Conductance Static Calibrations 

The properties of ordinary water are so complex that 
a theoretical conductance sensitivity is impractical. The 
important thing for these calibrations was not how do they 
compare with theory, but are the calibrations linear, 
repeatable, and measurable. 

Several wire conductance static calibrations were 
made using the static calibration technique discussed in 
chapter III-C-2 . One series of static measurements was made 
using the HP-4194A Impedance\Gain-Phase Analyzer with the 
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wires in either a ten gallon Pyrex beaker or the wind-wave 
tank. These measurements were made to compare the results 
from data taken in the two water samples. 

The conductance sensitivity of the measurements 
taken in the Pyrex is 6.4 ± 0.1 mS/cm. The conductance 
sensitivity of the measurements taken in the wind-wave tank 
is 6.3 ± 0.1 mS/cm. The data from both measurements are 
linear, as shown in Figure (4.3) . 
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Figure 4.3. Conductance vs. Relative immersion in two water 

systems (B3P31) . 

Another series of static wire measurements were made 
using the same equipment, in the wind-wave tank. These 
measurements were made to see if two' four-wire measurement 
combinations were possible using the five wire probe, and 
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how close the sensitivities were for the combinations. A 
gradient measurement technique was discussed, but not 
pursued. Table (4.3) shows the possible combinations. The 
wires that are used in each combination are shown filled in. 
TABLE 4.3. FOUR WIRE COMBINATIONS WITHIN THE FIVE WIRE PROBE 



(B3P39) . 



Four Wire Combination 


Four Wire Combinations 


Number 




SUl 




SU2 




SU3 


O##0# 


SU4 


•o### 



We used deviations in the conductance measurements 
to ensure that all the individual wires were straight as 
possible. We did use the tensioning loops to adjust some of 
the wires. The final measured conductance sensitivity and 
deviation results are shown in Figures (4.4 and 4.5) . 

The conductance sensitivities, shown in Table (4.4), 
of the various four-wire combinations are different because 
of the two different probe geometries, four in a row or a 
space inserted. These geometries put the voltage probes in 
different potential fields between the current probes. The 
SUl and SU3 configurations have lower sensitivities than the 
SU2 and SU4 configurations because the Av is the lowest in 
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this symetric configuration. The highest sensitivity should 
be obtained when the center probe is not used. 

The conductance sensitivity measurements were 
linear, with high correlation coefficients, and the 
deviations were consistent for all the combinations with a 
maximum error of 0.4%. The causes of the anomalies of the 
last data points for SU1 and SU4 are unknown. We choose to 
use wire combinations SU2 and SU4 in all latter test because 
they had the highest conductance sensitivities, as expected. 




R»Uiv* tmnwaion (cm). 



Figure 4.4. Conductance vs. Relative immersion for various 
four wire combinations (B3P42). 
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Figure 4.5. Deviation vs. Relative immersion for various 
four wire combinations (B2P42). 

TABLE 4.4. FOUR WIRE CONDUCTANCE SENSITIVITIES (B3P42). 



Four Wire 

Combination Number 


Conductance 

Sensitivity 

mS , „ „ . 
— (±0.1) 
cm 


SU1 


6.2 


SU2 


7.9 


SU3 


6.8 


SU4 


8.6 



Conductance system static calibrations were made to 
determine the voltage sensitivities which would become 
inputs for the engineering units of the HP-3562A Dynamic 
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Analyzer power spectrum measurements. The lock-in amplifiers 
provided electrical inputs to the conductance electronics. 
The conductance electronics electrical output was then 
measured using two hand held multimeters, one for each 
channel . 

The voltage sensitivity for channel one is 0.309 ± 
0.005 V/cm, and for channel two is 0.381 ± 0.005 V/cm. The 
data from the measurements was linear, as shown in Figure 
(4.6) . 




Figure 4.6. Output voltage vs. Relative immersion for 
conductance system electronics (B3P59). 
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B. DYNAMIC CALIBRATIONS 



1. Capacitance Dynamic Calibrations 

Several capacitance dynamic calibration apparatus 
and techniques were tried, as described in Appendix B. The 
apparatus was cumbersome, and the results of the techniques 
were not repeatable. The repeatability problem was caused by 
the low tank impedance, the results of this test were 
described in the previous chapter. 

It was also determined that all the capacitance 
dynamic calibration techniques were not conducted in the 
preferred motion either. These techniques moved the wire up 
and down in the water. This does not represent a traveling 
gravity wave motion. 

The low, approximately 20 ohm, tank impedance made a 
two wire capacitance system unusable when combined with the 
unknown and unstable water properties. Unfortunately, the 
tank impedance test was not conducted until late in the 
capacitance data collection efforts. 

2. Conductance Dynamic Calibrations 

Conductance dynamic calibrations were conducted to 
determine repeatability, if there was a preferred probe 
orientation, and to determine the frequency response. The 
calibrations were conducted in the small wave guide to cut 
off low frequency parametric waves, using the techniques 
described in chapter III-C-3. 
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Four consecutive measurements were taken in the same 
location, with three different probe orientations. The line 
of wires in the probe were either normal or parallel to the 
paddle. The first measurement of the sequence was always 
parallel to the paddle. Two measurements were always made 
with the wires parallel to the paddle. 

The two parallel measurements were the test of 
repeatability. The other measurements were the test for best 
probe orientation. The probe frequency response could be 
determined at all orientations. 

A two Hertz, one to four Hertz, dynamic power 
spectrum calibration was conducted in the small wave guide. 
Waves below two Hertz were not sinusoidal on the 
oscilloscope, which set the lower frequency limit. The Flat 
Top window was used for these measurements because this 
window has the lowest equipment, HP-3562A Dynamic Analyzer, 
amplitude accuracy uncertainty of up to ± 0.1 dB [Ref. 35]. 

The engineering units used in the HP-3562A internal 
power spectrum calculations were determined from the results 
of the static calibrations discussed in section A(2). 

The parallel probe orientation consistently had the 
highest sensitivity at any given frequency for both 
channels, as shown in Figure (4.7 and 4.8) . The data point 
at four Hertz is above the parametric wave cutoff for the 
small guide. The last data points show a definite 
interaction between the gravity and parametric waves for all 
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probe orientations. The parallel probe orientation was also 
very repeatable, except for the anomaly at the 2.5 Hertz 
knee . 

The HP-3562A Dynamic Analyzer measurements were 
verified by measuring wave peak-to-peak amplitudes using a 
HP-54501A Digitizing Oscilloscope. The amplitudes ranged 
from 330 mV at two Hertz to 187 mV at four Hertz (B3P68-90) . 

The very disappointing result was that the probe 
frequency response was not flat, not even for an octave. We 
have assumed that the probe response would continue to drop 
with frequency, but could not test the higher frequencies 
because of the parametric waves caused by the paddle. The 
current probe construction would not allow for a smaller 
wave guide, which would enable us to test at higher 
frequencies . 
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Figure 4.7. Sensitivity vs. Frequency for channel one 
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Figure 4.8. Sensitivity vs. Frequency for channel two 

(B3P64-67) . 

C . TANK MEASUREMENTS 

Measurements were taken in the wind-wave tank through- 
out the experiment. The sensor may or may not have been 
statically calibrated prior to any particular measurement. 
No tank measurements were made after the conductance system 
dynamic calibrations were made. 

Results of early phase speed and power spectrum 
capacitance system measurements were presented at the 119th 
meeting of the Acoustical Society of America, at Penn State 
University [Ref. 36]. 

The cross correlation plot in Figure [4.9] is a typical 
capacitance system measurement phase speed plot. The data 
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were analyzed using the least square fit technique. This 
1.35 ± 0.01 Hertz, driven gravity wave has a measured 
frequency of 1.36 ± 0.01 Hertz. The measured delay time is 
15.6 ± 0.1 msec. The probes had a measured separation of 
0 . 025 ± 0.01 m. 

The theoretical phase speed for this driven gravity 
wave, using equation [2.3], is 1.16 ± 0.01 m/sec. The 
measured phase speed, using equation [2.25], is 1.63 ± 0.01 
m/sec. This is a 40.7% difference. It is possible that 
electrical cross-talk between the two channels in the 
electronics box may have decreased the delay time measured. 

The theoretical wavelength for this same wave, using 
equation [2.26], is 0.86 ± 0.01 m. The measured wavelength 
is 1.20 ± 0.01 m. This is a 39.9% difference. 

The capacitance system equipment and techniques were 
both adjusted, and the measured phase speed and wavelength 
were consistently higher than theory. We believe the results 
are high due to an erroneous delay time in the channel two 
output caused by electrical cross-talk in the electronics 
box, but this was not confirmed. 

Base line noise and wind power spectrum measurement 
plots, using a statically calibrated capacitance system, are 
shown in Figure (4.10). Similar data obtained with the 
conductance system is shown in Figure (4.11). The four upper 
lines in Figure (4.10) are wind power spectrums, at four 
different wind speeds. Spectra for both the capacitance and 
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Hertz and 



conductance systems peak at approximately two 
decay at approximately 50 to 70 dB per decade, or with power 
decaying at a)~5 to or 7 , but the slope of the decay may not 

be trustworthy because the transducer was not dynamically 
calibrated and its frequency response is unknown. However, 
there is excellent agreement in the wind power spectrum 
measurements for the two transduction systems. 
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Figure 4.9. Capacitance system cross correlation plot . 
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Power spectra of (a) Base line noise, (b) driven gravity 
wave, and (c) driven gravity wave with wind, using a 
statically calibrated capacitance system, are shown in 
Figure (4.12). The driven gravity wave spectrum shows the 
harmonics of the driving frequency, at least 18 dB below the 
fundamental. There is very little energy above the noise 
level at other frequencies. 

The driven gravity wave and wind spectrum is 
approximately 20 dB above the noise and the driven gravity 
wave spectrum. The wind energy completely dominates the 
harmonics of the driven gravity wave with this particular 
frequency resolution. The wind decay is unaffected by the 
addition of the driven gravity wave. The decay of the power 
is approximately as 0 )" 6 . 

The power spectra plots of the driven gravity wave and 
driven gravity wave with wind were expanded to 0.5 Hz about 
the fundamental driven frequency, shown in Figure (4.13). 
The energy peak is downshifted 0.062 Hz when wind is added. 
A second, lower frequency energy peak is also observed. This 
peak is approximately 30 dB above the noise level. This peak 
may be caused by the non-linear effects discussed in Dr. 
Larraza's theory. Determining the cause (s) of the additional 
energy peak and the fundamental frequency downshifting were 
left for follow-on theses. 
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Figure 4.12. Capacitance system base line noise, driven 
gravity wave, and driven gravity wave with wind power 

spectra plot. 
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Figure 4.13. Expanded plot of driven gravity wave and driven 

gravity wave with wind. 
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Plots of the tank quality factor measured with the 
capacitance system, are shown in Figure (4.14) and similar 
data obtained with the conductance system is shown in Figure 
(4.15). The tank Q measurement with the conductance system 
was almost double the capacitance system measurement's value 
at the low frequencies, and half the value at the higher 
frequencies. It is unknown why and how there is a difference 
in tank Q for the two measurement techniques. The measured 
Q at different frequencies is shown in Figure (4.16). 

The driven gravity wave spectrum was not measured 
using the conductance system. All efforts were directed at 
getting a dynamic calibration for the conductance system 
instead. We were not sure of the best method to make phase 
speed and wavelength measurements with the conductance probe 
at that time. 
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Figure 4.14. Tank Q capacitance system measurement plot. 
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Figure 4.15. Tank Q conductance system measurement plot. 
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Figure 4.16. Tank quality factor vs. Frequency measured by 
capacitance and conductance systems (B2P137-142 and B3P10). 
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V. CONCLUSIONS 



Capacitance and conductance measurements using two and 
four wire techniques were developed, tested, and calibrated 
in this thesis. The four wire conductance system technique 
was determined to be the better of the two tested, and is 
currently being used in a follow on thesis. The two wire 
capacitance system repeatability problems caused by 
grounding and the complex properties of "simple" water are 
not present using the four wire technique. The two systems 
did give quite similar results for wind-wave spectra. 

We feel the four-wire technique is superior to the two- 
wire technique as long as the potential problems in 
conducting a water surface wave height measurement, dirty 
probes and unknown variable water chemistry, only exist on a 
microscopic scale. This thesis did not prove or set limits 
to this assumption. A study of these problems is suitable 
and desirable as a follow on thesis. 

Static calibrations were performed many times on all 
systems. The results of the various tests were consistent, 
linear, and generally compared well with theory. The static 
calibration techniques are simple and should be completed 
prior to other measurements in future experiments. 

The results of the conductance system dynamic 
calibration were disappointing and puzzling because the 
frequency response was not flat. It is difficult to infer 
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what the frequency response is in the entire frequency range 
of interest, 0-10 Hz, due to the limited results we were 
able to obtain. However, the two systems gave quite silimar 
results for the wind-wave spectra. 

Driven gravity wave frequency downshifting and wind 
energy dumping was observed in the combined gravity wave and 
wind-wave measurements. A follow-on thesis will explore 
this . 

The importance of the dynamic calibration is unknown for 
future nonlinear experiments and should be reviewed to 
insure that important data is available and interpretations 
are accurate. 

The most important thing I learned from this research 
are some basic principles of experimental science. First, 
you must have some idea concerning what and why you are 
investigating nature. Second, you must have a predetermined 
purpose and procedure for the investigation. Third, data 
must be recorded thoroughly, meticulously, and in a timely 
manner. Fourth, you must evaluate the data carefully to 
determine the outcome of your experiment. Finally, you need 
GOOD LUCK. 
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APPENDIX A 



I. HP-15C CALCULATOR LEAST SQUARES FIT STANDARD DEVIATIONS 

PROGRAM CODE 

The standard deviations in all the least squares fit 
parameters were discussed in chapter II E. This is program 
code for the HP-15C hand held calculator, which can be used 
to solve equations (2.16 through 2.19) for data entered into 
the calculator using the summation key (B1P66/67). The 
calculator stores the different summation values in the 
following register as shown in Table (A.l). 

TABLE A.l. HP-15C REGISTER VALUES AFTER USING SUMMATION KEY. 



Register 


Value 


2 


N, the number of samples 


3 


summation of x n 


4 


summation of x n ^ 


5 


summation of y n 


6 


summation of y 


7 


summation of x n Yn 



The standard deviations in y, slope, and intercept are 
displayed at program lines 040, 054, and 067 respectfully. 
The program code is shown below. 
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Program 
line number 



Calculator 

kevs 



Program 

function 

g[P/R] 

f [LBL] C 
RCL[5 ] 
g[x 2 ] 

RCL [ 4 ] 
x 
2 

ENTER 
RCL [7] 
x 

RCL [3] 
x 

RCL [ 5 ] 
x 

CHS 

RCL [2] 
ENTER 
RCL f 7 ] 
g[x 2 ] 
x 
+ 

RCL [2] 
ENTER 
RCL [4] 

X 

RCL [3 ] 
g [x 2 ] 



CHS 
RCL [ 6 ] 
+ 

RCL [ 2 ] 

ENTER 

2 



Vx 

STO [ 8 ] 
R/S 
RCL [4] 
ENTER 
RCL [4 ] 
ENTER 
RCL [2] 
x 



000 

001 


42,21,13 


002 


45,5 


003 


43,11 


004 


45,4 


005 


20 


006 


2 


007 


36 


008 


45,7 


009 


20 


010 


45,3 


Oil 


20 


012 


45,5 


013 


20 


014 


16 


015 


30 


016 


45,2 


017 


36 


018 


45,7 


019 


43,11 


020 


20 


021 


40 


022 


45.2 


023 


36 


024 


45,4 


025 


20 


026 


45,3 


027 


43,11 


028 


30 


029 


10 


030 


16 


031 


45,6 


032 


40 


033 


45,2 


034 


36 


035 


2 


036 


30 


037 


10 


038 


11 


039 


44,8 


040 


31 


041 


45,4 


042 


36 


043 


45,4 


044 


36 


045 


45,2 


046 


20 



107 



RCL [3 ] 
g [X 2 ] 



RCL [ 8 ] 
x 

R/S 
RCL [2] 
ENTER 
RCL [2] 
ENTER 
RCL [4] 
x 

RCL [3] 
g[x~2] 



RCL [8] 
x 

g [RTN] 



45,3 

43,11 

30 
10 
11 

45,8 

20 

31 

45,2 

36 

45.2 
36 

45,4 

20 

45.3 
43,11 
30 

10 

11 

45,8 

20 

43,32 



047 

048 

049 

050 

051 

052 

053 

054 

055 

056 

057 

058 

059 

060 

061 

062 

063 

064 

065 

066 

067 

068 
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APPENDIX B 



I. CAPACITANCE SYSTEM DYNAMIC CALIBRATIONS 

This section will briefly describe the equipment and 
techniques used for the capacitance system dynamic 
calibrations. Two dynamic calibrations of the capacitance 
system were explored before the capacitance system was 
determined to be unusable. Both of the techniques moved the 
water or the sensor probe, up and down, which we decided was 
not the preferred motion for the dynamic calibration. 

A . Cone Head 

A 16 ounce plastic funnel, or cone, was used to 
connect a 110 ml glass burette to a USRD Type J-ll 
Transducer, hence the Cone Head name for the apparatus. The 
J-ll provided a means to convert an electrical signal to 
fluid motion and the funnel acted as a fluid amplifier with 
a gain of approximately 130. All capacitance system 
measurements were actually made in the burette. 

The top of the probe wire was held in place at 
the top of the burette with a small cork stopper. Fishing 
line was tied to the bottom of the capacitance probe. The 
fishing line was passed around a large paper clip lodged in 
the funnel, pulled out of the top of the burette to tension 
the probe wire and held in place at the top of the burette 
with the same small cork stopper. 
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All joints in the Cone Head were sealed for 
water tight integrity with Polymeric Systems, Inc., PSI 601 
Clear, 100% Silicone, RTV Adhesive Sealant. The Cone Head 
apparatus is shown in Figure (B.l). 




Figure B.l. Cone Head apparatus. 

The J-ll transducer motion was calibrated using 
an LVDT. The J-ll compensation bladder pressure was measured 
using a PX-126, Omega Engineering, Inc. differential 
pressure transducer. This provided a non-intrusive 
measurement of the J-ll when the Cone Head was assembled. 

The electrical outputs of the LVDT and pressure 
transducer were compared when a low frequency square wave 
was applied to the J-ll. The Cone Head was assembled with 
water in the apparatus and the electrical output of the 
pressure transducer was then compared with the capacitance 
system output. 
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This calibration scheme was not refined due to 



the fact that only one sensor probe could be fit into the 
Cone Head at a time and the poor results in the Cone Head 
static calibrations. 

B. Motor 

This technique was briefly tested, with no 
documented results. The sensor probe was moved up and down 
in the wind-wave tank, while the probe was on the Harp. 

A 1/4 horsepower AC motor was bolted to the top 
U channel on the Harp. A 2.5 inch diameter pulley was 
mounted eccentrically to the motor shaft. A one inch 
diameter pulley was mounted directly above the probe. One 
end of an approximately 24 inch long, ten pound, fishing 
line was tied to the top of the probe, placed in the grooved 
flywheel, and the other end of the fishing line was tied to 
a standoff post on the 2.5 inch pulley. The other end of the 
probe was connected to the lower U channel with two number 
one rubber bands. 

The AC motor was plugged into a variable 
transformer which would adjust the mechanical output 
frequency of the large flywheel and provide the forcing 
function for the probe. 
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APPENDIX C 



I. HP-3562A DYNAMIC SIGNAL ANALYZER FRONT PANEL CONTROL 
SEQUENCE TO STORE DATA ON DISK 

The front panel hard disk accessing sequence for the HP- 
3562A Dynamic Signal Analyzer using the HP-9133 Disk Drive 
is shown in Table (C.l). Complete operating instructions can 
be found in Reference 35. 

TABLE C.l. HP-3562A HARD DISK ACCESSING SEQUENCE. 



Step 


Menu 


Select 


1 


HP-IP FCTN 


Select Address 


2 


SELECT ADDRESS 


Disc Address = 0 






Disc Unit = 1 


3 


DISC 


Save/Recall File 
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